par(pch=3, col="red") 
par(mfrow=c(1,1))
tenfile <- file.choose()
new <- read.csv(tenfile,header = TRUE)
n = nrow(new) 
n
#n = nchar(n)
atschool = 0
outOfSchool  = 0
for (i in 1:n)
{
  if(new[i,1] == "yes")
    atschool = atschool + 1
  else if(new[i,1] == "no")
    outOfSchool = outOfSchool + 1
}

#giả định số lần đi học tương ứng khả năng thành công
#vẽ biểu đồ xác xuất đi học của 1 VIên TRƯờng DHKHTN
p = atschool /(n-1)
k <- seq(0, (n-1), by = 1)
plot (k, dbinom(k, (n-1), p), type='l', ylab="Probability", col='blue',
      main = "Khả năng đến trường của 1 sinh viên DHKHTN",xlab="Số Ngày")